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We compare two strategies of multi-detector detection of compact binary inspiral signals, namely, 
the coincidence and the coherent for the realistic case of geographically separated detectors. The 
naive coincident strategy treats the detectors as if they are isolated - compares individual detector 
statistics with their respective thresholds while the coherent strategy combines the detector network 
data coherently to obtain a single detection statistic which is then compared with a single threshold. 
We also consider an enhanced coincidence strategy which is intermediate in the sense that though 
the individual statistics are added in quadrature and the sum compared with a single threshold, the 
estimated parameters are also checked for consistency. For simplicity, we consider detector pairs 
having the same power spectral density of noise, as that of initial LIGO and also assume the noise 
to be stationary and Gaussian. Further, since we consider the detectors to be widely separated on 
Earth, we take the instrumental noises to be uncorrelated; the wide separation implicitly means 
that since the detector arms must lie parallel to the Earth's surface, the detectors necessarily have 
different orientations. We compare the performances of the methods by plotting the receiver oper- 
ating characteristics (ROC) for the strategies. Several results are derived analytically in order to 
gain insight. Simulations are performed in order to plot the ROC curves. A single astrophysical 
source as well as a distribution of sources is considered. We find that the coherent strategy is supe- 
rior to the two coincident strategies that we consider. Remarkably, the detection probability of the 
coherent strategy is 50% better than the naive coincident strategy. One the other hand, difference 
in performance between the coherent strategy and enhanced coincident strategy is not very large. 
Even in this situation, it is not difficult to perform the real data analysis with the coherent strategy. 
The bottom line is that the coherent strategy is a good detection strategy. 

PACS numbers: 95.85.Sz,04.80.Nn,07.05.Kf,95.55.Ym 



I. INTRODUCTION 

Inspiraling binaries are one of the most promising 
candidates for first detection of gravitational waves 
(GW). The compact objects can be treated essen- 
tially as point particles leading to sufficiently ade- 
quate description of the system in terms of the post- 
Newtonian formalism [l[ . The great accuracy of the 
post Newtonian approximation of the phase, about 
a cycle for a wave train ~ 10 4 cycles long, renders it 
amenable for matched filtering analysis [2[ . Inspiral- 
ing binaries are astrophysically important, because 
they will not only carry detailed information about 
the binary system, but also general relativistic devi- 
ations from Newtonian gravity in their orbit can ex- 
perimentally be measured [1,3- The best available 
estimates suggest that at 1% false alarm probability 
the expected number of neutron star(NS)-NS binary 
coalescence seen per year by ground based interfer- 
ometers is 3 x 10~ 4 — 0.3 for initial detectors and 
1 — 800 for advanced detectors In recent years, a 
number of ground based detectors are taking qual- 
ity science data and are collaborating together, thus 
the time is ripe to consider analysis of network data 
for the detection of inspiraling binaries. The advan- 
tages of multi-detector search for the binary inspiral 
is that, not only does it improve the confidence of 



detection, but it also provides information about the 
direction and polarization state of the source. 

Two strategies currently exist in searching for in- 
spiraling binary sources with a network of detectors: 
the coherent and the coincident. The coherent strat- 
egy involves combining data from different detectors 
phase coherently, appropriately correcting for time- 
delays and polarization phases and obtaining a sin- 
gle statistic for the full network, that is optimized in 
the maximum likelihood sense. On the other hand, 
the coincident strategy matches the candidate event 
lists of individual detectors for consistency of the es- 
timated parameters of the GW signal. A coincidence 
search with real data has been carried out by several 
group [f| and data analysis on S5 data of the LIGO 
detectors has recently been performed Q- 

There is a long standing debate as to which strat- 
egy performs better. In an earlier paper 8], here- 
after referred to as paper-I, we compared the perfor- 
mances for the simple case of coaligned detectors 
located in the same place. The situation of two 
coaligned detectors located at the same place with 
correlated noise was first considered in paper-I and 
improved in is applicable in the case of existing 
detector pairs H1-H2. The coincident strategy has 
the advantage of reduction of false alarm, however 
this is at the price of reduced detection efficiency. 



2 



On the other hand, in the coherent strategy, the false 
alarm rate is not reduced, the sensitivity is enhanced 
which in turn results in higher detection efficiency. 
Which of these two competing effects wins was de- 
termined by looking at the detection efficiencies of 
the two strategies at the same false alarm rate. Re- 
ceiver operating characteristics (ROC) curve, which 
is the plot of detection efficiency versus false alarm 
rate, was drawn for both the strategies and from 
those curves it was inferred that for the viable false 
alarm regime the coherent strategy performs much 
better than the coincident strategy. 

In this paper we consider the general case of 
widely separated pairs of detectors. Since the detec- 
tors are situated on the globe and have their arms 
lying parallel to the Earth's surface, they necessar- 
ily have different orientations; we however consider 
the general problem of detectors with arbitrary lo- 
cations and orientations, since this does not greatly 
add to the mathematical complexity of the problem. 
The coherent statistic for non-aligned detectors is 
completely different from that of the aligned case 
[ill ]. We cannot extrapolate the results of paper-I 
to incorporate the more general case that we con- 
sider here. We compare the performance of the two 
strategies for the generalized case by obtaining the 
relevant ROC curves. We then investigate an en- 
hanced coincidence strategy which is basically an 
improvement on the simple coincidence strategy for 
it is seen that when the detectors are not aligned, 
the performance of the simple coincidence strategy 
is very poor. In enhanced coincidence, although the 
two detectors are considered in isolation and the can- 
didate event lists compared for consistency in the 
estimated signal parameters, the two statistics from 
individual detectors are added in quadrature and a 
single threshold is placed upon the resulting statis- 
tic. This strategy, for the case of two detectors only, 
renders a statistic identical in form to that of the 
coherent strategy. The details of this strategy are 
given in section III B.2. It should be noted that it 
is different from the coherent strategy. 

The plan of the paper is as follows. In section 
III1 we discuss the different coordinate systems re- 
quired to describe the problem, the signal and the 
response of the detectors and the network. In section 
IIII1 we discuss the false alarm and detection prob- 
abilities and derive the analytical formulae for the 
coincidence and coherent strategies. In section IIVI 
we perform simulations and then use these to plot 
the ROC curves for different detector pairs. We first 
briefly present the simulation method and then we 
determine the parameter windows required for coin- 
cidence analysis. We perform simulations for arbi- 
trarily oriented detector pairs and obtain the false 
alarm and detection probabilities for each pair and 
use these results to plot the ROC curves. In sec- 
tion |V1 we summarize our results and discuss future 
directions. 



II. THE SIGNAL AND THE RESPONSE 

In this section we discuss the GW signal and the 
response of individual detectors as well as of the net- 
work and introduce a normalization scheme and no- 
tations. We follow the conventions and notations 
of [ll|. In order to make the paper self contained 
we briefly review the formalism. It provides us with 
an efficient framework for the analysis here in the 
context of widely separated detectors with arbitrary 
orientations. 



A. Reference frames and Euler angles 

To understand the response to the GW signal of 
the detectors and the network, the first step is to 
identify the different coordinate systems naturally 
associated with the problem and the interrelations 
between them. The frames are related via rotations 
described by Euler angles which then appear in the 
network response in the form of beam pattern func- 
tions. 

Wave frame - (X, Y", Z) : The gravitational wave 
travels along the positive Z direction and X and Y 
denote the axes of polarization - (X, Y, Z) form a 
right-handed Cartesian coordinate system. 

Frame of detector I - (xx,yi,Zi) : This is the 
frame of the detector / having its origin at the inter- 
section of the arms of the detector. The arms of the 
detector lie in the xi — yi plane, which is the plane 
tangent to the surface of the Earth with the xi axis 
bisecting the angle between the two arms and yi is 
chosen such that the frame forms a right-handed co- 
ordinate system with the zi axis pointing radially 
out of the surface of the Earth. Note that the Earth 
is assumed to be a sphere for the purposes of this 
analysis. 

Earth frame - (x,y,z): This is the Earth frame 
with its origin at the centre of the Earth. The z axis 
points due north, the x axis points to the intersection 
of the equator and the Greenwich meridian and the 
y axis is chosen such that the coordinate system is 
right handed. This frame is usually chosen as the 
fiducial frame of reference, with respect to which 
the orientations and locations of each detector can 
be specified. In table U we specify the locations and 
orientations only of those detectors which we will use 
in our analysis. This is clearly not an exhaustive 
list of ground-based detectors existing or planned 
and the results obtained in this paper can also be 
obtained similarly for other detectors. 

Let (</>, 9, t/j) be the Euler angles through which 
the Earth frame must be rotated in order to align 
with the wave frame. The angles <fi and 9 are closely 
related to spherical polar coordinates and ip is the 
polarization angle. We opt for the Goldstein conven- 
tion [l2|. The Euler angles (ai , [3i that rotate 
the Earth frame to the frame of the J-th detector 
are then (A/ + tt/2, l u Cj + 3tt/4). 
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B. Signal at a detector 

Having described the coordinate system, we write 
down the response of the J-th detector to the GW 
signal in the Fourier domain. We adopt the following 
convention for the Fourier Transform, 



dt h(t)e 



Imft 



(2.1) 



as in paper-I. 

In the stationary phase approximation the spin- 
less, restricted post-Newtonian inspiral signal is 
given by, 

V{f) = N x £ J /- 7 /6 

xcxp j* 7 (/;t c ,<5 c ,r ,r 3 ), (2.2) 

where t c and 8 C are respectively the coalescence time 
and the coalescence phase of the binary. E 1 is the 
extended antenna pattern function which depends 
on seven angles - two angles (9, cj>) specifying the 
location of the source, two angles (i, ip) specifying 
the orientation of the source and three Euler angles 
{a i, Pi, 7/) specifying the orientation of the detector 
with respect to the wave frame. 



E 1 



1 + COS 2 I 



FL + i cos lFI 



(2.3) 



where the x are the usual antenna pattern func- 
tions depending on 9, 4>, tp, cij, 0i, 77 [Tl(. Note that 
E 1 is a complex quantity. Thus the phase of E 1 , say 
X 1 , where, tan% 7 = 2F^ cost/Fj_(l + cos 2 1) repre- 
sents the polarization phase and plays an important 
role in coherent detection. On the other hand we 
observe that x 1 simply adds on to the coalescence 
phase S c and therefore for the purposes of matched 
filtering in a single detector can be absorbed in 5 C . 
In the coincidence strategy each detector is treated 
separately, just like a single detector, and therefore, 
for matched filtering when using coincidence, x 1 can 
be absorbed in 5 C . 

The factor AT depends on the distance r to the 
binaries, the total mass M = mi + m-i and the re- 
duced mass ratio 77 — m ^J2 2 , where m\ and m 2 are 

Detector latitude(A) longitude(l') x-arm(£) y-arm 



H1,H2 46.45°N 

LI 30.56°N 

VI 43.63°N 

Kl 137.18°N 



-119.41°E 
-90.77°E 
10.5°E 
36.25°E 



36.8° 126.8° 

108.0° 198.0° 

71.5° 341.5° 

295.0° 25.0° 



TABLE I: The locations and orientations of some ground 
based GW detectors. LI - LIGO Louisiana, HI - LIGO 
Hanford 4 km, H2 - LIGO Hanford 2 km, VI - Virgo, 
and Kl - LCGT (the K stands for Kamioka). The lati- 
tude and longitude of the centre and the angles through 
which the x-arm and the y-arm must be rotated clock- 
wise (viewed from top) to point to North are listed. The 
angles (X,£, C) are related to (qj , f3i for a particular 
detector /. 



the individual masses of the binary. In the units of 
c = G = 1 it is given by, 

^^(Ay /2 ^ 5/ v /2 ^ 2/3 , M 

Instead of the total mass and the chirp mass of the 
binary system, it is customary to use the dynamical 
parameters tq and T3 for in these parameters the 
template placing is approximately uniform. These 
parameters are defined as: 



To 



256-Krifa 



(ttM/ q ) 



-5/3 



817/0 



(2.5) 



where f a is the fiducial frequency usually chosen 
to be the seismic cut-off frequency - the lowest fre- 
quency of the detection bandwidth. 

The phase of the signal 5' 7 (/; t c , 8 C , tq, T3) in the 
Ith detector frame relates to the phase in the fiducial 
frame *(/; t Cl 8 Cl t , t 3 ) as, 

*'(/; t c , 6 C , r , ts) = *(/; t c , 5 C , t , r 3 ) + 2nfAt I , 

(2.6) 

where At 1 is the time-delay between the J-th detec- 
tor and the fiducial detector. The fiducial detector 
may be taken to coincide with the origin of the Earth 
frame, that is, a detector situated at the centre of 
the Earth or it could be chosen as one of the detec- 
tors in the network. For the phase in the fiducial 
frame we adopt the same 3PN formula we took in 
paper-I (equations (2.7) and (2.8) in paper-I) given 

by 1- 



C. The matched filtering paradigm for a 
network of detectors 



By virtue of (|2.2p . we can represent the i-th. tem- 
plate for the J-th detector as, 

/i / (/;/2i,t c ,(S c ) = A 1 (so(/;Mi,*c) cos S c 

+ si /2 (f;fl t ,t c ) sin 5 C ) , (2.7) 



where, 



5 7r/2(/;/^A) = is (f;p h t c ) . (2. 



We require that both the templates Sq and s^ 2 have 
have unit norm; i.e. the scalar products (sq,Sq)i = 
( s i/2> s i/2) 1 = 1- The scalar product (a,b)i of two 
real functions a{t) and b(t) for the 7-th detector is 
defined as, 

(a,b) I = 2 df -j— , (2.9) 

where, we use the Hermitian property of Fourier 
transforms of real functions. Sl(f) is the one sided 
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power spectral density (PSD) of the noise. Although 
in this paper we do not consider the case of different 
PSDs for different detectors, we still keep the discus- 
sion general so that it could also be used for the gen- 
eral case of different PSDs. This does not complicate 
the discussion much, ft also follows from this defini- 
tion that when the templates are computed using the 
stationary approximation, (so,s^y 2 ) = ®> i- e -' they 
are orthonormal. Then, it is evident that A 1 is the 
amplitude of the waveform, i.e., (h 1 ,h T ) — (A 1 ) 2 , 
ft is the lower cut-off frequency which is normally 
taken to be the seismic cut-off - 40 Hz for initial 
L1GO - and is the the upper cut-off frequency 
usually taken to be about 1 kHz when the signal 
power for inspirals within the usual mass range, say, 
1M Q < m 1 ,m 2 < 40M© and initial L1GO PSD, 
falls of below a fraction of a percent. In most of the 
cases we have taken the noise PSD of the detectors 
to be identical and the index / can be dropped in 
Eq. (|2.9p . However, since for one case we consider 
different noise PSDs, we keep the subscript in this 
discussion to retain the generic nature. 
The statistic for the J-th detector is, 



the ratio of the total energy received by the system 
to the maximum possible energy accessible. 

This suggests the definition of the network- 
normalised signal to be 

S 1 = 7/6 exp [ **(/; t c , 6 C = 0, r 0j r 3 )] . 

9i 

(2.14) 

where, 



Q' = 



E 1 
\\E\ 



(2.15) 



The network vector Q = (Q\Q 2 , Q N ) lies in C N 
and has unit norm. This is the vector of polarization 
phases, which effectively 'brings all the detectors to 
the same orientation'. 

For the two detector network, the data consists 
of two data trains, {x*(t)\I = 1,2; and t E [0,T]} 
where data is taken in the time interval [0, T]. As- 
suming additive noise n 1 in each detector we have: 



1 = 1,2. 



(2.16) 



max 



(s I ,h I )j + (4/ a ,A J )? 



(2.10) 



which is to be compared with the threshold. 

Instead of two templates Sq and s^, 2 , we find it 
more convenient to use a single complex template 
S 1 which combines the two together: 

S x (f) = - r 7/6 exp [ »¥(/; t B , 6 C = 0, r , r 3 )] . 
9i 

The normalization factor gi is chosen such that 
(S^S 1 )! = 1. It is given by, 



9 2 i 



df 



P'*SUf) 



(2.12) 



In terms of <j/, the amplitude of the waveform h 1 is 
A 1 = ME 1 gi. 

In addition to the above normalization scheme fol- 
lowed for the construction of the template bank, par- 
ticularly for coherent detection, where the network 
is treated as a whole, we also need the concept of 
network normalization. The total energy accessible 
to a network of N detectors is a scalar and is given 

by 

N 2 



i=i 



7=1 



1-4 



NW 



(2.13) 



The quantity J2?=i 9i e *i eI = \ \ E \\ 2 is the C 2 norm 
of E 1 in C N . The quantity M 2 g 2 has the signifi- 
cance of the maximum possible energy accessible by 
the network i.e., the energy accessible when all the 
detectors are optimally oriented; so ||i?|| 2 denotes 



The noise random variables satisfy the statistical 
property: 

(n I (f)n}(f'))=ls I h (f)S(f-f), (2.17) 

where, the angular brackets denote ensemble aver- 
age. In this paper we take the noise in the two de- 
tectors to be uncorrelated, i.e, 



(« 1 (/K(/')> = o, 



(2.18) 



which is not unjustified if the detectors happen to 
be geographically separated. For the sake of sim- 
plicity of the analysis, we assume the noise to be 
both stationary and Gaussian. This is assumed in 
the analysis in section IIIII and for the simulations in 
section IIVI However, analysis of real data suggests 
it is neither - the SNRs and the \ 2 are correlated [? 
]• 

We conclude this section by introducing the defi- 
nition of the complex correlations C 1 which will be 
particularly needed for construction of the coherent 
statistic. We retain the notation and definition as 
in paper-I which was originally introduced in [rT| . 
The complex conjugate of C 1 is denoted by Cf. It 
is given by: 



CI = (£',: 



-7r/2 ■ 



(2.19) 



where, Cg and c^.y 2 are the real and imaginary parts 

of C 7 ; they are obtained by taking the scalar prod- 
ucts of the data x 1 with Sq and s 7 , 2 respectively; 
that is, 



4),7T/2 



( S 0,7t/2i X ) 



(2.20) 
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III. FALSE ALARM RATE AND 
DETECTION PROBABILITY 

Since in this paper we will be considering the 
case of two geographically separated detectors, we 
will consider pairs among the detectors HI, H2, LI, 
Virgo and LCGT, except for the pair HI and H2 
which happen to be in the same location. In this 
section we will derive and present analytical formu- 
lae for the false alarm rate and detection efficiency 
for the coherent and coincidence strategies. In co- 
incidence detection we consider two sub categories, 
one of straight forward coincidence which we call 
naive coincidence and a more sophisticated coinci- 
dence strategy which we call enhanced coincidence. 
We therefore organise our discussion as follows: 

A. Coherent detection 

B. Coincidence detection 

1. Naive coincidence detection 

2. Enhanced coincidence detection 

In the analytical formulae that we derive, we do 
not include all procedural steps of the numerical sim- 
ulation; instead, we take a simplified route. The de- 
rived formulae are thus slightly different from those 
obtained from the results of the simulation. Never- 
theless, the analytic formulae provide us guidelines 
by giving us functional forms for the false alarm and 
detection probabilities in terms of parameters which 
then can be determined from simulations. Even in 
such a simplified situation, as encountered in paper- 
I, where we deal with aligned detectors in the same 
location, the formulae depend on the unknown pa- 
rameters -/Vi n d - number of statistically independent 
templates and -number of statistically inde- 

pendent templates within the error window. These 
quantities are difficult, if not impossible, to evaluate 
analytically; we therefore determine them through 
simulations. We find here in the case of geographi- 
cally separated detectors, the quantity N w - U1 is signif- 
icantly different from that in paper-I, because here 
the time-delay window becomes larger by the maxi- 
mum time-delay between the two detectors. In this 
section, we obtain the relevant formulae in which the 
parameters N w i n and N{ n( i appear implicitly. The 
situation will be clear from the sections which fol- 
low. 



A. Coherent detection 

Coherent detection involves combining data 
streams in a phase coherent manner so as to ef- 
fectively construct a single, more sensitive detector. 
The maximum likelihood network statistic for two 
arbitrarily oriented geographically separated detec- 
tors has already been found. For the case of two 



arbitrarily oriented detectors the network statistic 
is given by [111 ]: 

A = ||C|| 2 = ic 1 ! 2 + |C 2 | 2 

= (clf + (ci /2 f + (^f + (cl /2 f, (3.1) 

where C 1 is the complex correlation of the I-th 
detector (7=1,2). 1 This is quite different from 
the coherent statistics for two coaligned detectors 
as it does not contain the terms involving cross- 
correlation between the two detectors. So, the false 
alarm rate and the detection efficiency changes con- 
siderably even for slightly non-aligned detectors. 
The comparison between the two strategies there- 
fore must be drawn separately for non-aligned de- 
tectors, for it is not only a more general situation, 
it is different altogether. In [l3| the two detector 
paradox regarding this abrupt change of detection 
statistic has been discussed and the improvement 
of the coherent strategy by suitably incorporating 
contributions from cross correlation terms has been 
advocated. We, however, restrict ourselves to the 
formalism of coherent detection as in [111 ]. Any im- 
provement to this basic formalism will enhance the 
performance of the coherent search and will further 
strengthen our results obtained in section [TV] 

Secondly, the correlations in the two detectors are 
computed at the same mass parameters {to, T3}, and 
at time-lags which differ at most by the light travel 
time d/c between the detectors, where d is the dis- 
tance between the geographically separated detec- 
tors. This is matched filtering in which the net- 
work template is matched to the network data. The 
statistic is the maximum taken over the permitted 
time-lags and polarization phases which make up the 
vector Q. We will see later that for enhanced coinci- 
dence, the same expression for the statistic appears, 
but the correlations in the two detectors are permit- 
ted to be evaluated at different points in the param- 
eter space, constrained by a certain window or error 
box (which will be specified later in the text). We 
therefore note that this is not matched filtering. 

To obtain the false alarm rate and the detection 

efficiency, first, we concentrate on a single template. 

1 2 

When there is no signal, Cq'^ , 2 are Gaussian random 
variables with zero mean and unit variance. So, the 
probability distribution of LLR can be shown to be 

Po(A) = jexp . (3.2) 



Note that in [TUl the square root of A is used. But here 
we use A because it has the x 2 statistics for Gaussian noise 
which has a simpler mathematical expression and thus eas- 
ier to implement. 



G 



The rate of false alarm for a threshold of A* is, 

dApo(A) 



pltcmplatc 

FA 



A* 



= (1 + T )exp 



A* 



. (3.3) 



We argued in papcr-I that when we have closely 
packed correlated templates, we can treat the tem- 
plate bank effectively as having Af;„d statistically in- 
dependent templates. iV; n< a is smaller than the actual 
number of total templates in the bank - it will turn 
out that A^i n d is substantially smaller. We will also 
make a simplifying assumption that the noise PSDs 
of the detectors are identical. With this assump- 
tion, the template placement is identical and so also 
is ATjnd. In which case, for the entire template bank, 
the false alarm rate is, 



Pfa = N ind 1 + 



A* 



exp 



A* 



(3.4) 



On the other hand, when, there is a signal of am- 
plitude -4nw in the network data, the expectation 
value of the squared statistic is A"^ w an( ^ the prob- 
ability distribution of the statistic is: 



Pi (A) = 



A 



1 



2 \ *4nw 



exp - A + ^ w MA N wVA). 



(3.5) 

where I\ is the modified Bessel function of the first 
kind. False dismissal occurs when in spite of the 
presence of the signal, the statistic A falls below 
the threshold A* . The false dismissal probability 
is given by, 



FD 



-f 

Jo 



dA ( VA 



r 2 \ Anw 
x h{A NW s/A). 



exp 



A 



A 2 



(3.6) 



The detection efficiency Pde or detection probabil- 
ity is then just P DE = 1 — P FD . 



B. Coincidence detection 

In coincidence detection, the two detectors are 
treated essentially in isolation. Separate lists of can- 
didate events are prepared; a candidate event occurs 
in a given detector I (where / = 1,2), when the 
statistic A/ = \C T \ 2 computed for the detector / 
crosses the threshold A} set for the detector /. This 
procedure produces two event lists, each for one de- 
tector. In absence of prior knowledge of the signal, 
we may choose the same threshold A* for both de- 
tectors. The next step involves matching the lists of 
candidate events and obtaining pairs of events. The 
matching is performed by ascertaining whether the 
estimated parameters for a pair of candidate events, 
each event chosen from a separate event list, lie in a 



predetermined parameter window. For a real astro- 
physical event, the estimated signal parameters must 
be consistent: ideally, the coalescence times t c must 
at most differ by the light travel time d/c, where d 
is the distance between the detectors, and the dy- 
namical parameters tq , t 3 must be identical. This 
ideal situation may only be realised in the limit of 
infinite SNR. However, for realistic SNRs, as are to 
be expected from astrophysical considerations, these 
constraints must be made less stringent, because the 
presence of noise introduces error in each parameter. 
The estimated parameters therefore will differ from 
their actual values; they must lie within an error- 
window W. The determination of the size of this 
window is pivotal to this strategy. 

The window W is determined as follows. Consider 
first the parameter t c . For two geographically sepa- 
rated detectors, separated by distance d, the maxi- 
mum time-delay is d/c. The estimated values of t c in 
the two detectors can differ in addition to d/c from 
errors due to noise. We denote the error boxes in 
t c due to noise by A„t c (Ai) and A„i c (A 2 ) in detec- 
tors 1 and 2 respectively. The error boxes depend on 
Ai and A 2 for the respective events and also on the 
probability we can tolerate in losing an event. We 
choose the probability of not losing the event to be 
99 % for each of the three parameters, giving a final 
probability of not losing an event to be (.99) 3 ~ 0.97 
for the error window. The error box A„f c (Ai, A2) is 
determined numerically by carrying out simulations 
(see next section) . We assume that the errors in the 
two detectors are independent and so the total error 
Ayyic is then realised as a quadratic sum. We have: 



(A w t r y = J-j + (A„t c (Ai)) 2 + (A„t c (A 2 ))^. 

(3.7) 

Similar considerations hold for the dynamical pa- 
rameters r , t 3 and their corresponding error boxes: 

(A w r ) 2 = (A„T (A 1 )) 2 + (A„t (A 2 )) 2 , (3.8) 

and 

(A w r 3 ) 2 = (A„r 3 (A 1 )) 2 + (A„r 3 (A 2 )) 2 . (3.9) 

Let At c , Ato, At 3 be the differences in the estimated 
parameters for the candidate pair of events corre- 
sponding to the parameters i c ,To,T 3 respectively, 
then we say that the parameters match if, 

I At c I < Awt c , I At I < A w t , |At 3 | < A w t 3 . 

(3.10) 

Note that the dimensions of the window W arc 
symmetric under the interchange of the detector la- 
bcllings 1 and 2. This is important for consistency: 
if a trigger in detector 1 has a window which includes 
a trigger in detector 2, then the trigger in detector 
2 has a window of the same size, albeit translated, 
which then includes trigger 1. Let Af w i n (Ai,A2) be 
the number of templates in the window now deter- 
mined by both Ai and A 2 . Note that because of the 
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symmetry in definition of the window size, we have 
AT wln (Ai,A 2 ) = iV w i n (A 2 ,Ai). 

Recently, a geometrical method for determining 
W using Fisher information matrix has been pro- 
posed which uses ellipsoidal windows instead of a 
rectangular ones [l^]. This scheme takes into con- 
sideration the correlation between the parameters 
which then reduces the false alarm rate. Incorpo- 
rating ellipsoidal windows would have the effect of 
shifting the ROC curve for coincidence strategy a 
little to the left. We may estimate by how much the 
curve would shift. From the figure 1 shown in 0] of 
the ellipse in the (to, T3) plane describing a template 
one can see that the area of the ellipse is about a 
quarter of the corresponding rectangle. However, for 
the parameter t c the window must take into account 
the time-delay for geographically separated detec- 
tors. This is done by the translating the error el- 
lipsoid (due to noise only) around a trigger in one 
detector by the amount ±d/c and check whether the 
so translated ellipsoid intersects the error ellipsoid of 
the second detector. If d/c is reasonably larger than 
the errors in t c only due to noise, then there is al- 
most no advantage to be gained in this parameter in 
using ellipsoidal windows. Thus the false alarm may 
be reduced by about a factor of 5 which is ~ 0.7 
on the logarithmic scale (the logarithm is to base 
10). Although the procedure would produce better 
results, for simplicity of implementation we perform 
the conventional coincidence search with rectangu- 
lar windows and introduce a fudge factor a ~ 1/5 
in 7V W i n which takes into account the effect of using 
ellipsoidal windows in reducing the false alarm. 

1. Naive coincidence 

The procedure for implementing naive coincidence 
for two arbitrarily located and oriented detectors is 
as follows: 

1. Choose the same threshold A* for the two de- 
tectors. Let A/, I = 1,2 be the individual 
statistics of the two detectors, prepare two can- 
didate event lists such that A/ > A*, 7=1,2. 

2. Look for pairs of candidate events, each can- 
didate event coming from a different list, such 
that the sets of estimated parameters match, 
that is, the condition (|3.10|) is satisfied for the 
two events. 

If the above requirements are satisfied, announce 
detection. We will now use these conditions to ob- 
tain expressions for the false alarm rate and false 
dismissal. 

In paper I we have shown that for a network of two 
detectors, the expression for the probability of false 
alarm was obtained assuming a fixed window size 
dependent only on the threshold A*. Here we have 
a variable window size now depending on the SNRs 



pi (or cquivalcntly A/) of the two events and hence 
the derivation is more involved. However, we nor- 
mally need the result at high value of the threshold 
A* in which case the expression for the false alarm 
probability assumes a simple form. 

Let the statistics Ai,A2 > A* cross the thresh- 
old A*. The probability that A 2 lies in the in- 
finitismal interval of size <iA 2 for some given tem- 
plate is ^e~ A2 ^ 2 dA2- Now this is the approxima- 
tion we make: we assume that there is exactly one 
false alarm in detector 2 lying within the window; 
in general we could have more than one false alarm 
in the window, but if we assume a high value of the 
threshold A*, then it is unlikely that one has more 
than one false alarm occurring within the window 
and thus this probability can be neglected (usually, 
in the literature the probability of at least one false 
alarm is calculated, but finally it is usually approxi- 
mated to the probability of exactly one false alarm; 
here instead we directly compute this probability). 
Thus in this approximation, the probability of A 2 
lying in an interval [A 2 , A 2 + dA 2 ] for the templates 
in the window is N win (A u A 2 ) x ±e~ A2/2 dA 2 . This 
probability must be multiplied by the probability 
A^ind x ^e _Al / 2 (iAi to obtain the probability in the 
rectangle [Ai, Ai + dAi] x [A 2 , A 2 + dA 2 ] and finally 
integrated from A* to 00. Thus for the probability 
density we obtain the expression: 

p(Ai,A a ) = iiV ind iV win (A 1 ,A 2 )e-( Al + A2 )/ 2 . 

(3.11) 

The false alarm probability Pfa is then given by in- 
tegrating this probability density over the acceptable 
region. Thus, 

/>oo />oo 

P FA = / dAi / dA a p(Ai,A 2 ). (3.12) 
Ja* Ja* 

For high threshold A* and also using the fact that 
Awi n (Ai, A 2 ) falls off rapidly with increasing Ai, A 2 , 
it is easily seen that this expression approximates to, 

P FA ~ N ind N win (A*, A*)exp (-A*) . (3.13) 

As remarked before, iV w i n depends crucially on 
the threshold A*, because the size of the error win- 
dow strongly depends upon the SNR as suggested 
by both Fisher information matrix considerations as 
well as by the simulations we perform. It also de- 
pends weakly on the location of the signal in the 
parameter space. However, we choose to ignore this 
weak dependence and treat N W i n as a function of 
the threshold only as in paper-I. Exact value of N w - m 
must lie between zero and the total number of tem- 
plates within the parameter window. We choose 
N w - m to be the total number of templates within the 
parameter window following the practice adopted in 
paper-I and such a choice leads to good agreement 
between theoretical estimates and numerical simula- 
tions. 
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We now turn to the detection probability. When 
a signal of network amplitude Anw falls on the de- 
tectors, the expected value of the statistic in each 
detector will be different because of differing orien- 
tations of the detectors. This information is encoded 
in the extended antenna pattern functions E 1 which 
depend on the detector orientations. Since the noise 
in the detectors is uncorrelated, the detection effi- 
ciency is just the product of the detection efficiencies 
of the individual detectors: 

ft E =PrWPrk (3.14) 

where, Pl, E , 1—1,2 are the detection efficiencies of 
the individual detectors. In earlier literature, these 
have been calculated and here we merely state the 
result: 

^de = i£°dAcx P j-i[A + (.4 / ) 2 ]j 

x IoiA 1 ^!). (3.15) 

Note that each A . , J = 1,2 is proportional to the 
extended antenna pattern functions E 1 respectively. 
The presence of the individual functions E 1 in the 
expression of detection efficiency demonstrates the 
weak point in this type of coincidence detection. In 
general, the response of the two detectors will be 
different because of different orientations, but since 
one has no a priori knowledge of the direction of 
the source, one sets the same threshold for the two 
detectors which leads to a lot of false dismissal. In 
fact, for the case when the configuration is such that 
the response of one detector is close to zero for some 
particular direction to the source, actual candidate 
events would be missed even if the noise level is low. 

The main reason behind the poor performance of 
the naive coincidence strategy has its roots in the sky 
coverage. For a particular source at a fixed distance, 
the intrinsic SNR of the GW emitted depends cru- 
cially on the location and orientation of the source 
(the angles 9, 0, -0 and t). For a single detector, and 
for a given threshold (minimum SNR) , such a source 
would be "visible" for certain directions and orien- 
tations. When the detectors are aligned the sky cov- 
erage of the coincident detector is the same as any 
one of the single detectors although the coincident 
detector still performs better than the single detec- 
tor only because of low false alarm rate. For the 
non-aligned case, since different parts of the sky are 
covered by each detector, for naive coincidence, it is 
the intersection of the sky coverage of the two indi- 
vidual detectors which decides detection. This cru- 
cially hampers the naive coincidence strategy. Just 
to get a quantitative idea of the performance of this 
strategy, we consider a binary source at a distance 
15 Mpc. The mass of each individual star is 1.4M Q . 
The binary is taken to be optimally oriented, i.e., 
t/> = i = 0. The location of the source in the sky 
is varied and the intrinsic SNR is calculated for LI, 
VIRGO orientations and LIGO I design sensitivity. 



In a simplistic scheme, if the intrinsic SNR in a given 
detector is above 7 for a particular direction, we take 
that direction of the sky as covered. We find that the 
coverage of a single detector is 49% while in naive 
coincidence the detectors cover only 18% of the sky. 
The corresponding number for the coherent strategy 
is 92% of sky coverage. 

In order to counter the disadvantage of low sky 
coverage of this simplistic strategy, we consider a 
more sophisticated coincidence strategy which over- 
comes this problem. We describe this strategy in the 
next subsection. 



2. Enhanced coincidence detection 

As we have pointed out for non-aligned detectors, 
the performance of coincidence strategy is poor be- 
cause of low sky coverage. We adopt a more sophis- 
ticated strategy which overcomes this problem. The 
procedure is as follows: 

1 . Choose a low threshold Ag and if A/ , 7=1,2 
are the individual statistics of the two detec- 
tors, prepare two candidate event lists such 
that A/ > Ao\ 1= 1,2. 

2. Look for a pair of candidate events, the events 
coming from separate lists, such that the sets 
of estimated parameters match within the 
error-window; the window has already been 
defined in Eq. (|3~T0| . 

3. Choose the final (high) threshold A* > 2A^ 
and construct the final statistic A = Ai + A2 
and register detection if A > A* . 

This strategy in essence leads to increased detec- 
tion efficiency. Further, this strategy also involves 
preparing separate candidate lists and matching can- 
didate events, but for the case of two detectors, the 
detection statistic is of the same form as that of the 
coherent strategy. The detection statistic is, 

A = ||C|| 2 = \C 1 \ 2 + \C 2 \ 2 

= (4) 2 + (4/ 2 ) 2 + (^) 2 + (4 /2 ) 2 . (3.16) 

Note that although this statistic has the same form 
as in the coherent case, as remarked before, it is ac- 
tually different: the mass parameters here do not 
have to be the same for the two detectors, but 
are only constrained to lie in the window described 
above. Also in t c , the error box is somewhat larger 
than in the coherent case, where the time-lag must 
differ at most by d/c; in addition to d/c there is er- 
ror introduced by noise which increases the size of 
the error box. 

We first obtain the false alarm probability. We 
first note that the first two steps in the procedure are 
identical to those of naive coincidence except with 
the final threshold A* for naive coincidence replaced 
by the low threshold Aq applied in the first step. 
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Therefore, identical derivation follows, and hence the 
probability density here is the same as before and 
is given by Eq. (|3.11[) . However, the integration 
region is different essentially because of the third 
step. Define the regions (or events) Ri,R 2 ,R\ 2 in 
the (Ai, A2) plane: 

Ri = {(A 1 ,A 2 )|A 1 > A*, 0<A 2 <c»}, 
R2 = {(Ai, A 2 )|A 2 > Aq, 0<Ai<oo}, 
R l2 = {(Ai,A 2 )|Ai + A 2 > A*}. 

(3.17) 

Let the region R = R\ n R 2 D R\ 2 be the intersection 
of all the three regions. Note that since A* is cho- 
sen greater than 2Aq we have R as a proper subset 
of Rx fl i?2- The region 'cut out' from i?i n R 2 is 
an isosceles right angled triangle, we call A. Thus 
Ri H R2 = R + A where + denotes disjoint union of 



the sets. The false alarm probability is given by in- 
tegrating the probability density given by Eq. (|3.1ip 
over the region R. Thus: 

Pfa- / dA 1 dA 2P {A 1 ,A 2 ). (3.18) 

JR 

We may be able to approximate this expression by 
writing R = R\ n R 2 — A which leads to: 

P FA = P(Ri n R 2 ) - P{A) . (3.19) 

If again N w - m is a rapidly decreasing function of both 
its arguments, we can again write: 

P(Rx n R 2 ) w N ind N win (A*, A*)e~ A '° . (3.20) 

On the other hand, P(A) becomes 



Thus, we have 



P(A) 



= / dx\dx 2 


p(xi,x 2 ) 




J A 






N ind r A '- 
4 4 


-AS 

dx\ e 






Ja* 


~ N ™« N 

— 4 wm 




/•A* -AS 

/ dx 

K 



dx 2 N win (x 1 ,x 2 ) e 



A^ind AT win (A*,AS) 



e 2 



A* 



Ao + 1 



(3.21) 



P FA ~ N ind N win (A* , A^ e" 



A* 



-Ao + 1 



(3.22) 



We now turn to the detection probability. We 
assume a signal in the detectors with amplitudes 
A 1 , 1—1,2 and true signal parameters Aq , where 
A p = {t c , tq, T3}. The true signal parameters {tq, T3} 
will be the same for the two detectors while the pa- 
rameter t c will differ by the light travel time depend- 
ing on the direction of the wave. Thus Aq must be 
interpreted according to the context. However, be- 
cause of noise in the detectors, the detector statistics 
Ai and A 2 in general will be evaluated at points in 
the parameter space different from Aq , say at, A^' and 
Ag , which will lie close to the true signal parameters 
Aq, if the signal amplitude is sufficiently high. We 
now compute the probability density function (pdf) 
for the statistic Ai in the limit of high amplitude 
for detector 1. Similar arguments will hold for the 
pdf for detector 2. In the limit of high amplitude 
we may assume the errors in the estimated parame- 
ters in detector 1 to be Gaussian distributed with a 
covariance matrix Cf. Also the amplitude param- 



eter of the signal will be be reduced by the ambi- 
guity function Tl(AX^), where we assume that the 
ambiguity function only depends on the difference 



in signal parameters AA^ 



An. The ambi- 



guity function is normalised to TL(0) = 1 which is 
also its maximum value. (Since we have assumed 
identical noise PSD for the two detectors, the am- 
biguity functions for the two detectors are the same 
and accordingly we have omitted the detector index 
I for the ambiguity function.) The joint pdf for the 
statistic Ai and signal parameters A± is given by, 



p 1 (X^A 1 ) = f(X^g(X^A 1 ), 



(3.23) 



where, 



(27r) 3 / 2 (detCi)5 



■ i[cr 1 UAA5'AA;' 



(3.24) 
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and, 



ff(Af;A x ) 



ig-£[Ai+(.A 1 ) 2 W 2 (AA<»)] 
2 



IoiA 1 W(AA?) v^T]. (3.25) 



A similar expression obtains for detector 2 with the 
index 1 replaced by the index 2. 

We now note that the covariance matrix scales as 
the inverse of the square of the amplitude or Cj* ~ 

CI" (A 1 = I) /{A 1 ) 2 and in the limit of A 1 > oo, 

the function /(Af) tends to a delta function centered 
at the signal parameters Aq . Now we integrate over 
the signal parameters A^ or marginalise over them. 
Assuming high amplitude A 1 , we expand g up to the 
second order in AA^ around Ag, and the result after 
integration is: 



Pi(Ai) =pi(AA^ =0;Ai) + 



fc^CTC^ = 1) 
(A*) 



(3.26) 

where fc* are constants, namely, the second deriva- 
tives of g with respect to AA M evaluated at AA M = 0. 
This means that the pdf for detector 1 is essen- 
tially given by the first term plus corrections of order 
o(l/(A 1 ) 2 ). These corrections are small if A 1 >> 1. 
Or writing this result explicitly: 

Pi (AO = ie-^+^l/o^v/Ai) +o(l/(A 1 ) 2 ) . 

(3.27) 

For the detector 2 we obtain an identical result with 
1 replaced by 2 on the RHS. Finally, the pdf for 
the two detectors is obtained by simply multiplying 
the pdfs of the individual detectors since the statis- 
tics Ai and A2 are independent. Thus the detection 
probability or efficiency is obtained by integrating 
the product of the pdfs over the region R and is 
given by, 

Pde= /*pi(Ai)pi(A a ) dAxdA 2 . (3.28) 

JR 



In this strategy the sky coverage is better than the 
naive coincidence because AJ is chosen much smaller 
than the threshold A* chosen in naive coincidence 
strategy. Thus we expect the performance to be bet- 
ter than the naive coincidence. In the next section 
we compare the performances of all the strategies. 



IV. COMPARING PERFORMANCES 



In the previous section we have derived analyti- 
cal formulae for plotting the ROC curves. However, 
the quantities, N- m & and N w - m , are unknown and we 
need to determine them by numerical simulation be- 
fore we can actually plot the curves. On the other 
hand, as we shall see below, the false alarm probabil- 
ity and the detection efficiency from the numerical 
simulation do not necessary agree with theoretical 
formulae perfectly due to various practical reasons. 
Thus, we have to resort to numerical simulations for 
evaluating quantitatively the performances of each 
detection strategy. Nevertheless, the analytical for- 
mulae provide guidance towards deriving useful fit- 
ting formulae from the simulation which are then 
used to plot the ROC curves. In paper-I we have 
already discussed this matter in detail. However, 
we include a description for making the paper self 
contained. 



Simulation method 



In this paper, we take the noise PSDs of all the de- 
tectors to be identical. We use the design sensitivity 
of initial LIGO [lj] for all the detectors. 



S h (f) = S [(4.49a;)- 56 + 0.16a;- 
= 00 otherwise, 



0.52 + 0.32a; 2 



/ > Is 



(4.1) 



r 



where, x = f/f k , fk = 150 Hz, f s = 40 Hz, and S = 
9 x 10~ 46 Hz -1 . For computing the SNR (scalar 
product) we set fi = f s . We take the upper cut-off 
frequency /„ to be 1024 Hz and the data are sampled 
at 2048 Hz in accordance with the Nyquist theorem. 
This choice of the noise PSD is simply to save the 
computation time of the simulation. If we take the 
PSD of advanced LIGO, the lower cut off frequency 
f s becomes around 10 Hz. The length of the signal 
in the bandwidth becomes longer which makes the 



computation time of the matched filtering longer. 
Even if we were to adopt the PSD of advanced LIGO, 
we do not expect much qualitative change in the 
results. We consider 4 detector pairs listed in Table 

El 

Using 3PN restricted waveform, we place the tem- 
plates 25 ms apart in the To direction and 19 ms 
apart in the T3 direction. This corresponds to a max- 
imum mismatch of 3%. We have taken a rectangle in 
the parameter space with 25 templates along each of 
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Pair Light travel time [ms] 



Ll-Vl 


26.42 


Ll-Hl 


10.02 


Ll-Kl 


32.47 


Vl-Kl 


29.18 



TABLE II: The light travel time between pairs of detec- 
tors. 



the coincident conditions, namely, Eq. (|3.10[) and sec- 
ondly A 2 (j) > Aq. If such a A 2 (j) exists, the pair 
{A]^), A 2 (j)} define the trigger list of the naive coin- 
cidence detection. From this trigger list, we compute 
a enhanced coincident statistic defined as 

A = A 1(i) +A 2W . (4.6) 



the To,T3 directions, totally giving us 625 templates 
in the ro,r3 parameters. We take a data train 32 
seconds long and sample it at 2048 Hz. We generate 
32 seconds Gaussian data and repeat the coherent 
and coincident searches for -/V s i m number of times. 
We test for various values of N s i m , and we find that 
N a [ m = 20000 is sufficient to give reliable estimates 
of the false alarm rate. Therefore the total data 
length of the simulation is 6.4 x 10 5 seconds. 

In the coherent searches, we first compute the 
statistics from each detector, Ai(£ c , r , t 3 ) = (cj) 2 + 
( C i /2 ) 2 and A 2 (t c ,T ,T 3 ) = {elf + (cl /2 ) 2 (i.e., 
Eq. (|3.ip ). which are functions of t c ,To and T3. We 
compute Ai + A 2 and take the maximum over the 
possible time delay and over the mass template 
space. We obtain the following statistic, 

A(t c ) = 

max max [Ai(i c , To, T3) + A 2 (t c + td, To, T3)] , 

T0,T3 \t d \<d/c 

(4.2) 

where d is the distance between the detectors and 
c the speed of light. Since the coalescence time t c 
is sampled at 2048 Hz, the samples A(t c ) (sampled 
at the same rate as t c ) are correlated. In order to 
remove the correlation, following the procedure in 
Paper I, we divide the data train into St = 15.6ms 
short data trains, where the i— th short train starts 
at time tu\ — 15.6ms x i. In each short data train, 
we take the maximum of A(£ c ), 

A (i) = max A(t c ). (4.3) 

<*c<*(i+l) 

{A^)} defines the trigger list of the simulation. 

In the coincident analysis, we first compute the 
statistics of each detector and take the maximum 
over the mass template space. As in the coherent 
analysis, we divide the data train into St — 15.6ms 
short data trains, and take the maximum of the 
statistics in each data train. We obtain the triggers 
from each detector defined as 

A 7(i) = max A/(t c ), (4.4) 

t(i)<t c <t(i+V) 

A/(i c ) = maxA 7 (t c ,r ,r 3 ), (/ = 1, 2) . (4.5) 

To,r 3 

In the naive coincidence detection, for each A^,) > 
Aq, we select a maximum of A 2(j -) which satisfies 



B. Coincidence windows 

To determine the size of windows for coincident 
analysis, one way is to use the Fisher information 
matrix for this purpose. However, it has been shown 
that at low SNR, SNR < 10, the Fisher information 
matr ix g rossly underestimates the size of the win- 
dow [l5j. So, we determine empirically the number 
of templates falling within the error- window W such 
that the signal is detected 99% of the time. Since we 
have three parameters, the coalescence time t c and 
the dynamical parameters tq and T3, on the average, 
the fraction of signals detected will be (.99) 3 ~ 0.97. 
The simulations are carried out for a single detec- 
tor for SNR(= \/~Kj) ranging from 4 to 25. Deter- 
mination of W at low SNR is difficult because the 
signal tends to get overwhelmed by false alarm trig- 
gers - in most cases, instead of the signal we pick 
up false alarms. Since the false alarm triggers oc- 
cur anywhere in the parameter space at random, it 
becomes impossible to decide on the window. So, 
we restrict ourselves only to those signals for which 
SNR > 4. For SNR> 25, we take the same window 
size as SNR= 25, since the parameter estimation ac- 
curacy becomes nearly constant at large SNR due 
to the finite mesh size of the mass template space 
and the finite sampling rate. This constant win- 
dow size at SNR> 25 does not affect our analysis 
very much because we will use the false alarm prob- 
ability and the detection probability corresponding 
A(= SNR 2 ) < 200. The parameter estimation errors 
from the simulation are summarized in Table Hill It 
is useful to deduce a fitting formula of Table IIII1 We 
observe that the estimated error scales roughly in- 
versely as the SNR which we denote by p. Thus, we 
may write: 

A„i c = % (4.7) 
VA 

A„r = ^, (4.8) 

A„r 3 = % (4.9) 
VA 

for p < 25 with a tc = 98.7ms, a TQ — 1895.4ms, 
and a T , = 1386.0ms. The window size is then 
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CMP 


A + 


A n T() 




4 


24.2 


4/0 


Q 1 /I 

314 


5 


23.1 


413 


289 


6 


19.6 


351 


260 


7 


16.0 


300 


234 


o 
8 


11.0 


2.39 


188 


9 


9.8 


194 


157 


10 


8.0 


162 


129 


11 


O.O 


lot) 


lUo 


12 


5.9 


126 


101 


13 


5.3 


110 


89 


14 


4.9 


109 


87 


15 


4.8 


98 


80 


20 


4.1 


82 


70 


>25 


3.6 


79 


66 



TABLE III: Average errors in parameters t c , to and T3 
in units of milliseconds due to detector noise for SNR= 
vAj ranging from 4 to 25. The size of the error box 
corresponds to 99% detection. For SNR> 25, we take 
the same window size as SNR= 25. 



given assuming the box window as: 



( A w t c ) {A w t ) ( A w r 3 ) 



c 2 ' 
Ai 



Ai 

< 
A 2 



a t 
A 2 
n l/2 r 



1/2 



,1/2 



T3 

Ai 



T 3 

A 2 



,(4.10) 



where d is the distance between the two detectors. 
As discussed in paper-I, the decorrelation length of 
t c , To and T3 can be taken as ~ 15.8ms, 220ms and 
80ms. The number of independent templates in the 
coincidence window size is thus estimated roughly as 



Nl 



(0) 


JLf. 


win 


at \ 


" 1 


1 


_Ai A 2 



1 

aT 



1 

T2 



1/2 



N, 



(o) 



15.8ms 220ms 80ms 



932.4. 



(4.11) 



(4.12) 



Note however that we do not use the value of 
N^ n explicitly. We use only the functional form of 
Eq. (|4.11[) to obtain the formulae for the false alarm 
rate. 



C. Estimating the false alarm rate 

The number of false alarm triggers as a function 
of the threshold from simulations is shown in Fig[TJ 
In these figures, and in the rest of the figures in this 
paper, "coincident (AND)" and "coincident (SQS)" 
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FIG. 1: The number of false alarm triggers as a function 
of the threshold for various detector combinations. 
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indicate the naive coincident strategy and the en- 
hanced coincident strategy respectively. We extrap- 
olate these plots to larger threshold in order to ob- 
tain the false alarm probability corresponding to a 
much larger data set. In this extrapolation, the ana- 
lytic formulae in the previous section are helpful. In 
the single detector case, the number of false alarms, 
]Vfa(A > A*) is given analytically as, 



iV FA (A > A*) = N ind e 



A* 



(4.13) 



The factor in the exponential, namely, —1/2 does 
not necessary hold in a simulation. We thus replace 
— 1/2 and also lniV in( j with parameters a and b to 
be determined from the simulations. We take the 
logarithm of the above equation and write: 



]nN FA (A > A*) = -aA* + b , 



(4.14) 



where a is the slope and b is the intercept if we plot 
lniVpA(A > A*) versus A*. From simulations we 
plot the required curve and read off a and b. 

We apply similar methods to coherent and coinci- 
dent cases. In the analytical formula for the coherent 
case, the number of false alarms falls off as Eq. (|3.4[) : 



In jV FA (A > A*) = -^A* + In (\ + ^ j + lnN ind , 

(4.15) 

In the case of naive coincidence, the number of false 
alarms is given by (Eq. (|3.13p h We use Eq. (|4.1ip for 
iV w i n and obtain, 



ln7V FA (A > A*) = -A* + - In 



1 (d 

c 



2 

A* 



+ ln(^J+ln(AWV^). 

(4.16) 

In the case of enhanced coincidence, we have from 
Ec, rf3~m 



In A^fa(A > A*) 



1 /A* 

-2 A * + Ht- a 5 + 1 

+ ln(N ind N win (A*,A*)) , 

(4.17) 

where Aq = 4 2 . The dominant term which deter- 
mines the slope of each curve is — A*/2 in Eq. (14.15)) , 
-A* in Eq.gTlD, and -A*/2 in Eq.gTT]). We 
replace these factors with a, and unknown terms 
lnM nd in Eq.glSl), ln(N iad N^l) in Eq.g3U) and 
ln(A^i n( jiV w i n (Ao, A )) with b respectively. We then 
determine a and b from the simulations. By fitting 
the curves obtained from simulations we determine 
a and b in the four cases: (i) the single detector 
case, (ii) the coherent case, (iii) the naive coinci- 
dence case and finally (iv) the enhanced coincidence 
case. The results are listed in Table II VI in this or- 
der for the four detector pairs (LI, VI), (LI, HI), 



Ll-Vl 


Fitting region 


a 


b 


Single 


30 < A* < 50 


-0.44567 


23.5293 


Coherent 


40 < A* < 60 


-0.4551 


26.6683 


Naive coincidence 


23 < A* < 33 


-0.89681 


34.0749 


Enhanced coincidence 


50 < A* < 67 


-0.46412 


31.2254 


(a) Ll-Vl 


Ll-Hl 


Fitting region 


a 


b 


Single 


30 < A* < 50 


-0.44624 


23.5448 


Coherent 


40 < A* < 60 


-0.4668 


26.3345 


Naive coincidence 


23 < A* < 33 


-0.84920 


33.1193 


Enhanced coincidence 


50 < A* < 67 


-0.46104 


31.0549 


(b) Ll-Hl 


Ll-Kl 


Fitting region 


a 


b 


Single 


30 < A* < 50 


-0.44596 


23.5328 


Coherent 


40 < A* < 60 


-0.45420 


26.8005 


Naive coincidence 


23 < A* < 33 


-0.90159 


34.0837 


Enhanced coincidence 


50 < A* < 67 


-0.46407 


31.2229 


(c) Ll-Kl 


Vl-Kl 


Fitting region 


a 


b 


Single 


30 < A* < 50 


-0.43849 


23.2738 


Coherent 


40 < A* < 60 


-0.45480 


26.7335 


Naive coincidence 


23 < A* < 33 


-0.90073 


34.1248 


Enhanced coincidence 


50 < A* < 67 


-0.46410 


31.2301 


(d) Vl-Kl 



TABLE IV: The results of the fitting of the number of 
false alarm triggers in various detector pairs. 



(LI, Kl) and (VI, Kl). We can see that in all cases, 
the slope a is slightly smaller than theoretical value. 
This is mainly because of the maximization over the 
coalescence time and the two mass parameters de- 
scribed in Eq. (|4.2[) - (|4.5[) . When we take the maxi- 
mum of A, we tend to pick up large A events. This 
produces the deviation from the simple theoretical 
curve. We have confirmed this by performing simu- 
lations in which no maximization is done. When no 
maximization is done, we obtain the slope a which 
agrees with the theoretical value in all of cases from 
(i) to (iv). In this paper, the maximization about 
t c , the time delay, and two masses are essential to 
eliminate the statistical correlation between nearby 
templates. Thus, the deviation of the slope a from 
the theoretical value is not unnatural. Further, in 
the analysis with real data, this type of maximiza- 
tion process is usually carried out. Thus, such a 
deviation is more realistic than the values obtained 
from simple theoretical analysis. 



D. The ROC curves 

After ascertaining the formulae for the false alarm, 
we proceed to plot the ROC curves. We first con- 
sider a hypothetical source such that the mass of 
each star in the binary is equal to 1.4M Q . We con- 
sider the detector pair (LI, VI). We take the source 
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to be located at 35 Mpc with the inclination angle 
and the polarization angle as 0. The direction to the 
source is taken to be described by the polar angles 
9 = 1.0 rad and <fi = 0.8 rad in the Earth centered 
coordinate system. The amplitude of the signal in 
the noise free situation is then given by A 2 = 32.2 for 
LI and A 2 = 28.3 for VI. The detection probability 
for this case is shown in Fig[2J Note that the de- 
tection probabilities in Fig[2]are slightly worse than 
the theoretical values obtained in Section III. This 
is because, in the simulation, we use a discrete time 
step and also discrete set of mass parameters values 
which then produce a mismatch between the signal 
and templates. Thus, the amplitude of the signal 
detected is smaller than injected value. 

From the detection probability and the false alarm 
rate obtained in the previous subsection, we plot 
the ROC curves. We use the fitting formulae for 
the false alarm rate. We assume an year's worth 
of observation period and templates in each individ- 
ual mass ranging from 1M Q to 40Af Q . The total 
number of time samples is then 6.1 x 10 10 and the 
number of mass templates is 1.2 x 10 4 . On the other 
hand, the number of mass templates in the simula- 
tion is 625 and the total number of time samples is 
32 x 2048 x 2 x 10 4 = 1.3 x 10 9 . Thus, the false alarm 
probability for one year observation period with tem- 
plate mass range of 1 — 40Af Q is found by scaling the 
false alarm probability by the factor: 



(a) Ll-Vl 

Detection Probability (a fixed source case) 



1.2 x 10 4 6.1 x 10 1 



625 



1.3 x 10 9 



= 900.9. 



We therefore add ln(900.9) = 6.8 to b in Table [TV] 
to the relevant Ll-Vl part of the Table. As dis- 
cussed in Section III.B, the coincidence window for 
two coincident strategies can be improved using el- 
lipsoidal windows with a factor of 5 reduction in the 
false alarm rate. We thus subtract ln(5) from b for 
naive and enhanced coincidence strategies in Table 
W\ We show the ROC curves in Fig. [3] We find that 
the coherent case gives the largest detection proba- 
bility. It is about 70 to 80% larger than the naive 
coincident case. The relative difference between the 
coherent case and the enhanced coincident case is 
about 10 to 20%. 

In order to obtain a fair comparison of the strate- 
gies, we consider a distribution of sources with differ- 
ent position, orientation and distance - we compare 
their average performance. The sources are taken to 
be uniformly distributed within 15 Mpc. The orien- 
tation and direction of the binary is randomly chosen 
from an uniform distribution (uniform in ip, cos t, (j) 
and cosf?). The detection probability is shown in 

Fig. a 

We again consider a one year data train and mass 
ranges from 1M Q to 40M Q . The ROC curves are 
drawn for Ll-Vl, Ll-Hl, Ll-Kl and Vl-Kl network 
in FigO A reduction of the factor of 5 for coinci- 
dent strategies has been included in these plots. For 
the sake of comparison, single detector performance 
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Single(L1) 
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150 



FIG. 2: The detection probability for coherent, coinci- 
dence and single detector for a fixed source in the LIGO 
(LI) and VIRGO (VI) case. The source is a 1.4-1.4 M Q 
binary located at 35 Mpc with ip = 0, l = 0, 9 = 1.0 rad 
and 4> = 0.8 rad. The amplitude of the signal at each 
detector is A = 32.2(L1) and A = 28.3(V1) respectively. 
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FIG. 3: The ROC curves for coherent, coincidence and 
single detector for a fixed source in the LIGO (LI) and 
VIRGO (VI) case. The source is a 1.4-1.4 M Q binary 
located at 35 Mpc with ip = 0, i = 0, 6 = 1.0 rad and 
<p = 0.8 rad. The amplitude of the signal at each detector 
is A = 32.2(L1) and A = 28.3(V1) respectively. 



curves are also drawn. We have assumed that all the 
detectors have the noise PSD of initial LIGO. 

For example, in the Ll-Vl case, at the false alarm 
probability of 10~ 5 , the detection probability is 0.69 
(coherent), 0.45 (naive coincidence), 0.66 (enhanced 
coincidence), and 0.45 (single). Thus, the coherent 
strategy performs around 50% better than the naive 
coincident strategy. In paper-I for coaligned detec- 
tors at the same place, we concluded that the co- 
herent detection is better than the coincident detec- 
tion by about 25-40%. Thus the difference in mis- 
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FIG. 4: The detection probability for coherent, coinci- 
dence and single detector for distributed sources. Two 
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FIG. 5: The ROC curves for coherent, naive coincidence 
and enhanced coincidence for (a) Ll-Vl, (b) Ll-Hl, (c) 
Ll-Kl and (d) Vl-Kl for distributions of sources. 
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aligned detectors is larger than the coaligned case. 
In the case of misaligned detectors, the amplitude 
of the signal in each detector differs, because of the 
different orientations of the detectors. This differ- 
ence produces smaller detection probability in the 
naive coincident strategy. This is manifest in Fig. [4j 
Moreover, the detection probability of the naive co- 
incident strategy is even less than each of the single 
detector cases of LI and VI. In the ROC plot (Fig. 
[5]), however, the detection probability of the naive 
coincident strategy and that of the single detectors 
is nearly equal because the false alarm probability of 
the naive coincident strategy is smaller than the sin- 
gle detectors which compensates for the difference 
in the detection probability. The reason why the 
naive coincidence's curve and single detectors' curve 
in Fig. [5] coincide is only by chance or coincidence. 
In the Ll-Hl case, for example, the light travel time 
is smaller and the orientations of the detectors is 
nearly the same and therefore, the naive coincident 
strategy gives larger detection probability than the 
Ll-Vl case for a given false alarm probability. 

Enhanced coincidence strategy performs better 
than the naive coincidence strategy but we see that 
coherent strategy is still superior by around 5%. Al- 
though the detection probabilities of the coherent 
and enhanced coincident strategies are nearly equal 
as can be seen in Fig. QJ the false alarm probabil- 
ities of the two strategies differ. This produces a 
difference in the ROC curves. 

Note that although in the two coincident strate- 
gies, we reject the injected signal with Aj < 16, 1 = 
1, 2 in the evaluation of the detection efficiency, we 
do not follow this procedure for the coherent strat- 
egy and in the single detector cases. Thus, if we take 
a threshold less than 16 in naive coincidence and 32 
in enhanced coincidence, we cannot do a fair com- 
parison. However, in the ROC curve of FigO the 
threshold range used in plotting the ROC curve is 
35 — 70 (naive coincidence) and more than 100 (en- 
hanced coincidence). Thus, the rejection of the low 
amplitude signals in the evaluation of the detection 
efficiency does not affect the ROC curves. 



V. CONCLUSION 

We compare the two strategies for analysis of net- 
work GW data for inspiraling binaries, namely the 
coherent and the coincident strategy. Analysis with 
real data will lead to actual comparison between the 
two strategies but as suggested by our findings under 
simplifying assumptions, the performance of the co- 
herent strategy is superior to the coincidence strat- 
egy on the whole. We distinguish two sub strategies 
in coincident detection, namely, naive coincidence 
and enhanced coincidence. For naive coincidence, 
the difference in performance, as compared with co- 
herent, is even more conspicuous for misaligned de- 
tectors than coaligned detectors, for there are re- 



gions of the sky from where the signal may not be 
detected separately by the two detectors ruling out 
the possibility of coincidence detection, while, with a 
coherent detector, the source may still be visible. In 
fact this difference is so glaring that another strategy 
which we call enhanced coincidence needs to be de- 
vised. The coherent strategy uses the statistic which 
is optimal in the maximum likelihood sense for the 
network. It inherently incorporates the information 
about the phase to decide on the detection for the 
statistic explicitly contains data from different de- 
tectors added with consistent phases. On the other 
hand coincident strategy treats the detectors sep- 
arately, missing the crucial information about the 
phase altogether in consequence. This vital differ- 
ence leads to superior performance of the coherent 
strategy. 

Although the coherent strategy is superior to co- 
incident strategies, the difference between the co- 
herent and enhanced coincident strategies is not so 
large. Only a relative improvement of about 5% in 
the detection probability is obtained with the co- 
herent strategy. One may ask whether there is any 
practical advantage in using the coherent strategy in 
the case of two misaligned detectors. Note however 
that the coherent method is not so computationally 
expensive compared with two coincident methods, 
since we do not take cross correlation of two detec- 
tors' data in the coherent strategy. Thus overall, the 
coherent strategy is a good detection method. 

In future we would like to consider detector net- 
works with more than two detectors. In case of co- 
herent strategy, adding a detector to the network 
increases the sensitivity and therefore the detection 
efficiency of the coherent detector without exception. 
The locations in the sky for which the source is de- 
tectable also increases for coherent strategy upon 
addition of detectors thereby increasing the number 
of potential sources. In fact regardless of sensitivity 
of a detector in a network, adding a detector to the 
network, for coherent strategy, always improves the 
performance of the network. Despite the increase in 
computational cost, conducting coherent search will 
enhance the probability of detection of GW. The sit- 
uation for enhanced coincidence is not so clear and 
needs investigation. 
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